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Abstract. It is shown that all contracting, spatially homogeneous, orthogonal 
Bianchi cosmologies that are sourced by an ultra-stiff fluid with an arbitrary and, 
in general, varying equation of state asymptote to the spatially flat and isotropic 
universe in the neighbourhood of the big crunch singularity. This result is employed to 
investigate the asymptotic dynamics of a collapsing Bianchi type IX universe sourced 
by a scalar field rolling down a steep, negative exponential potential. A toroidally 
compactified version of M*-theory that leads to such a potential is discussed and it is 
shown that the isotropic attractor solution for a collapsing Bianchi type IX universe 
is supersymmetric when interpreted in an eleven-dimensional context. 
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1. Introduction 

The possibility that our universe underwent a 'pre-big bang' contraction before bouncing 
into its present expansionary phase continues to attract attention (see, e.g., [HI21E1)- 
In scenarios of this type, a central question to address is the behaviour of the 
universe during the final stages of the collapse in the vicinity of the big crunch 
singularity. Recently, Erickson et al. jl] have considered the class of collapsing, spatially 
homogeneous cosmologies that are sourced by an ultra-stiff perfect fluid with an equation 
of state p = (7 — where 7 > 2. In the case of a constant equation of state 
parameter, they have proved that the universe asymptotes to the spatially flat and 
isotropic Friedmann- Robert son- Walker (FRW) cosmology on the approach to the big 
crunch. This is to be expected on qualitative grounds, since the energy density of the 
matter source scales as p oc a" 37 , whereas the curvature and anisotropics vary as or 2 
and a -6 , respectively. The latter therefore become subdominant as the spatial volume 
a? — > 0. More generally, for a variable equation of state, they have argued that if the 
spatial curvature and anisotropy are initially small, they remain so during the collapse. 

In this paper we provide a rigorous proof of a 'cosmic no hair' result for all 
contracting, orthogonal, spatially homogeneous Bianchi cosmologies, where the matter 
source has an arbitrary equation of state 7 = j(p) that is subject only to the condition 
that it is differentiable and bounded such that (7 — 2) is positive-definite for all time 
(including at the big crunch). One example of a matter source with an equation of 
state 7 > 2 is a minimally coupled scalar field with a negative-definite self-interaction 
potential. We employ the cosmic no hair theorem to gain insight into the behaviour of 
a collapsing Bianchi type IX cosmology sourced by a scalar field interacting through a 
steep, negative exponential potential. We find that such a universe can indeed isotropize 
at the big crunch. 

Negative exponential potentials are known to arise in compactifications of higher- 
dimensional theories such as M*-theory j^J. This is an eleven-dimensional theory in a 
spacetime with signature (9+2) and is directly related to M-theory [Hj by string dualities. 
For a particular toroidal compactification of this theory to four dimensions, we interpret 
the isotropic attractor solution for a contracting Bianchi type IX cosmology in an eleven- 
dimensional context and find that it corresponds to a supersymmetric background. 

2. Cosmic No hair Theorem 

2.1. Bianchi Cosmology 

Bianchi models are spatially homogeneous cosmologies admitting a three-parameter 
local group G 3 of isometries that acts simply transitively on spacelike hypersurfaces S t . 
Coordinates can be chosen such that the four-dimensional line element has the form 
ds 2 = —dt 2 + h a b(t)uj a u! b (a,b = 1,2,3), where the one-forms uo a satisfy the Maurer- 
Cartan equation du a = \C a \ }( pj h A uj c and C a b c are the structure constants of the Lie 
algebra of G 3 . Since C a ^ c ) = 0, C a b c has at most nine independent components and 
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these are classified in terms of a symmetric 3x3 matrix, n ab , and the components of a 
3x1 vector A b = C a ab . This implies that the structure constants can be expressed in 
the form 

C c ab = n cd e dab + 5f a A b] , (1) 

where e abc is the totally antisymmetric tensor with e 123 = 1. Substitution of Eq. ([Tj) 
into the Jacobi identity C e d[ a C d bc ] — then implies that A b is transverse to n ab : 

n ab A b = 0. (2) 

If A b ^ 0, it represents an eigenvector of n ab with zero eigenvalue and it may be assumed 
without loss of generality that n ab = diagfni, n 2 , n 3 ] and that A b = (A, 0,0). Moreover, 
a suitable rescaling can always be found such that the eigenvalues of n ab take values 
{0, ±1}. As a result, Eq. (J2J) simplifies to 

7ii>4 = 0. (3) 

Bianchi class A models satisfy A = and the class B are defined by the property A ^ 
(ni = 0) 0. 

Throughout this paper, we consider orthogonal Bianchi models where the fluid 
velocity vector is orthogonal to the group orbits. In the case of a scalar field matter 
source, this implies that the field is constant on the surfaces of homogeneity. We employ 
the orthonormal frame approach developed by Wainwright and collaborators [HIGH]. (For 
a review see In this approach, the Hubble scalar is defined by H = £/£, where I 

is a length scale and a dot denotes differentiation with respect to t. The Raychaudhuri 
and Friedmann equations are then given by 

H = -H 2 - \a 2 - 1(37 " 2)p (4) 
3H 2 = a 2 -^R + p, (5) 

respectively, where a 2 = \a ab a ab is defined in terms of the shear tensor a ab , is 
the scalar curvature of the t = constant hypersurfaces, and p and p represent the 
energy density and pressure of the fluid. The Friedmann constraint equation (JHJ) may 
be expressed in the form 

VL + T? + K =1 (6) 

by defining expansion-normalized density, shear and curvature parameters: 

p a 2 {3) R 

Q ^3IP> ^ = SIP' K= ~6IP- (7) 
A deceleration parameter, q = —££/£ 2 , may also be defined such that 

H = -(l + q)H 2 (8) 

and Eqs. (jlj) and (JHJ) together imply that 

g = 2S 2 + ^(3 7 -2)a (9) 
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Finally, the covariant conservation of energy-momentum is expressed in the form 
of the fluid equation: 

p = -SHjp, (10) 

where p = [j(p) — l]p defines the equation of state parameter, ■y(p). We will consider 
an arbitrary equation of state subject only to the conditions that p = p(p) is at least C l 
and that (7 — 2) is positive-definite for all values of the spatial volume. We also assume 
the standard energy conditions hold and, in particular, that p > 0. 

The scalar curvature satisfies < for all Bianchi types except the type IX. 
It then follows from the Friedmann equation (jSJ) that the spatial volume of an initially 
contracting Bianchi type I- VIII universe will decrease monotonically with cosmic time, 
t. In view of this, it proves convenient to define a dimensionless time variable r |%j ll()j: 

- = -- (11) 

Hence, r is a monotonically decreasing function of t when H < and, since < £ < 00, 
the big crunch singularity will occur at r — > — 00. 

We now show that Q —>■ 1 in the vicinity of the big crunch for all Bianchi type 
I- VIII universes. Firstly, Eq. (JBJ) implies that Q < 1 since K > 0. Furthermore, Eqs. 
©? ©) © an d ([TO"]) yield an evolution equation for the density parameter which takes 
the form: 

Q! = [-(3 7 - 2)K + 3(2 - 7 )S 2 ] VI, (12) 

where a prime denotes d/dr. It follows from Eq. (|T2")) that Q' < for any initially 
contracting Bianchi type I- VIII universe, with equality iff K = S 2 = for any non- 
vacuum orbit (Jl > 0). Hence, Q is a monotonic decreasing function of r and, since Q 
is bounded, we may conclude that lmv^-oo Q' = 0. Eq. (|T2|) then implies that 

lim K = lim S = 0, lim 0=1, (13) 

r — > — 00 r — > — 00 r — > — 00 

where the latter limit follows directly from Eq. (JHJ). 

To proceed further, we will require the specific form of the Einstein field equations 
for each Bianchi model. For the Bianchi type I- VIII universes, these equations can be 
expressed in the form of an autonomous set of ordinary differential equations (ODEs) 
of the form x' = f (x), subject to a constraint equation g(x) = 0, where the state vector 
xgK 6 . The physical interpretation of the state variables x is different for the class A 
and B models, however, and we therefore consider each class in turn in the following 
Subsections. 



2.2. Bianchi Class A (I - VIII) 

The physical state of a Bianchi class A cosmology is determined by the vector x = 
(H, S + , £_, Nx, N2, A/3), where N a = n a /H, E± = a±/H and a± are linear combinations 
of the two independent components of the (traceless) shear tensor j^J IBj- Thus, £± 
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determine the anisotropy associated with the Hubble flow and are related to the shear 
parameter by 



(14) 



whereas N a parametrize the spatial curvature of the group orbits and are related to the 
curvature parameter such that 



K = — 
12 



(15) 



Nf + Ni + N£ - 2(N 1 N 2 + N 2 N 3 + N 3 N X ) 

The Einstein field equations for the Bianchi class A are then given by an autonomous 
set of first-order ODEs jH] : 

= -(2 - q)Z± - S± (16) 

(g-4E + )7V 1 (17) 

(g + 2S + + 2v / 3S_)iV 2 (18) 

(g + 2S + -2v / 3S_)AT 3 , (19) 



where 



S + 



1 



2y/3 



(iv 2 -iV3) 2 -iv 1 (2iv 1 -iv 2 -iV3; 

(Ns-NJiNt-Nz-Ns), 



(20) 
(21) 



together with the decoupled equation 

H' = -(l + q)H. (22) 

Eqs. (11 H|l and (fT4^1 imply immediately that rim r _ > _ 00 S + = lim T ^_ 00 XL = and it 
follows from Eq. Q that lim T __ OQ q = (37 — 2)/2 > 0. By following a similar argument 
to that of ^2] (who consider expanding cosmologies with < 7 < 2/3), we may then 
deduce from Eqs. (fT7j) - (fT§j) that for each N a there exists a parameter e > such that 
N' a /N a > e for a sufficiently negative r, and hence that lim T ^_ 00 N a = 0. We conclude, 
therefore, that the spatially flat and isotropic FRW universe is the global sink for all 
ultra-stiff, orthogonal, initially contracting Bianchi models of type I- VIII (class A). 



2.3. Bianchi Class B 

For the class B models, we employ the framework developed by Hewitt and Wainwright 
[TO] and define the quantities a = \<j a ba ab and a + = \o a a , where a a b is the trace-free 
part of the shear tensor a a b- Likewise, we define h = ^h a bh ab and n + = \n a a) where 
n a b is the trace- free part of n a t. Since A 7^ for this class, there exists a constant h 
such that h = \{n 2 + — hA 2 ). For n 2 n 3 7^ 0, this defines the group parameter h = h^ 1 
for types VI^ {h < 0) and VII^ {h > 0), respectively (Types IV and V correspond to 
h = and the Bianchi type III is the same as the type VI_i). 

The specific form of the Einstein field equations is unimportant in establishing 
the nature of the attractor solution for ultra-stiff, collapsing Bianchi class B universes. 
Indeed, it is sufficient to define Hubble normalized variables |H1 HD] such that S + = 
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(J+/H, £ = a/H,A= A 2 /H 2 , N + = n + /H and N = §(jV* - hA), where A > 0, E > 0, 
and A > 0. In this case, the evolution equation for the density parameter, Q, still has 
the form given by Eq. (|T2|) . where the deceleration parameter, q, is defined by Eq. (JHJ), 
but the variables {K, E 2 } in these expressions are now defined by 

K = N + A (23) 

E 2 = Y? + + E. (24) 

Inspection of Eq. ()13|) therefore implies that as the universe approaches the big crunch, 
liniT-^-oo E + = lim T ^_ 00 E = and, since iV and A are both non-negative quantities, it 
also follows that lmv^-oo iV = lmv-^oo A = 0. Moreover, the dimensionless curvature 
variables for the Bianchi class B are given by 

S + = 2N, S ab S ab = 24(i + Nl)N, (25) 

from which we deduce that lim T ^_ 00 S + = limT^.oo S ab S a b = 0. 

In the Bianchi class B, there is a special model corresponding to the type VI_i/ 9 , 
which is exceptional in the sense that its isometry group does not necessarily admit an 
Abelian subgroup that acts orthogonally transitively. Nonetheless, Eqs. (|T2|) and ffTT^f) 
and Eqs. (|2^ j) - (j23j) remain valid for the exceptional type VI__i/ 9 and the above analysis 
therefore applies to all Bianchi class B models. (More precisely, the shear parameter 
(|2^|l acquires additional terms due to the extra independent component, but it may still 
be expressed as a sum of non-negative terms that must all vanish when Eq. (|T3*|) holds 
[13]). We may conclude, therefore, that the spatially flat and isotropic FRW universe 
is the global sink for all ultra-stiff, orthogonal, initially contracting Bianchi class B 
cosmologies. 



2.4- Bianchi Type IX 

It now only remains to consider the Bianchi type IX universe. The physical state of a 
Bianchi type IX model is parametrized by the state vector x = (H, a + , o~-,ni, n 2 , n 3 ) G 
3ft 6 as for other types in the class A. However, since n x > 0, n 2 > 0, and n 3 > for this 
model, the curvature parameter K is no longer semi-positive definite and, consequently, 
the Hubble parameter may not necessarily be a monotonically varying function of time. 
This implies that alternatives to the Hubble-normalized variables are required for a 
global analysis. A set of appropriate variables that compactifies the Bianchi type IX 
phase space was introduced by Hewitt, Uggla and Wainwright (see section 8.5.2 of [TT]). 
In this approach, physical quantities are normalized in terms of a function, D, such that 

(tf, £± ,^),(§,-- JL), (26) 

where 



D = ^H 2 + ^( ni n 2 n z fl\ (27) 
In particular, Eq. (J27|) implies that H is bounded by the constraint equation 

H 2 + i(A!A 2 A 3 ) 2/3 = 1. (28) 
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A new time variable, r, is also defined: 

dt 1 , . 

Tf = D < 29) 

and Eqs. (fTTjl - lfiTjj) . (J22j) and (J27j) then imply that the evolution equation for D takes 
the form 

D* = -(l + q)HD, (30) 

where a star denotes differentiation with respect to f and q = H 2 q. It can be further 
shown that the Einstein field equations take the form ^T] 

H* = -{1-H 2 )q (31) 

E* ± = -(2 - q)HE± - S± (32) 

N* = (Hq - 4E+)Ni (33) 

N; = (Hq + 2S+ + 2V3E_)iV2 (34) 

AT* = (Hq + 2S+ - 2V3£_) N 3 , (35) 

where S± are defined by Eqs. (j2Uj) and (|2*T|) . respectively, with N a replaced by N a . 
Moreover, it follows from Eqs. (jHJ) and (j2Hj) that 

1 -(3 7 -2)(l-y) + ^(2- 7 )S 2 , (36) 



where 



and 



2 v , , 2 



si + s 2 (37) 



V = — 
12 



iV{ + A 2 + ATg - 2NiN 2 - 2N 2 N 3 - 2jVijV 3 

+3(A 1 A 2 A 3 ) 2/3 ] . (38) 

The definition (J38)l implies that V > and substitution of Eq. (|T5j) yields V" = 
H 2 [K + (N1N2N3) 2 / 3 /4}. Thus, the Friedmann constraint equation © may be expressed 
in the form: 

S 2 + y + = l, (39) 
whereas substitution of Eq. (|3U|) into Eq. (JBTIj) implies that 

g = 2S 2 + ^(37-2)a (40) 

Finally, we may derive an evolution equation for the density parameter, Q = 
p/(3D 2 ), by differentiating with respect to f, and substituting in the fluid equation 
(fTUjl and the evolution equation (|3*D|). We find that 

n* = £lH [-(37-2)^ + 3(2 -7)S 2 ] . (41) 

We may now deduce from Eq. (|4()jl that when (7 — 2) is positive-definite, q > 
with equality iff fl = T? = 0. Furthermore, since Eq. (|2*Hj) implies that H is bounded, 
— 1 < H < 1, it follows from Eq. (|3*Tj) that H* < and, consequently, that 5" is 
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a monotone decreasing function. On the other hand, q > for a non-vacuum orbit 
(Q > 0), and this implies that \im f _^ 00 H — 1 and lim f ^ +00 H = — 1. In other words, an 
initially expanding model will eventually undergo a recollapse (when H passes through 
zero). Once the recollapse occurs, Eq. (}4~T|) implies that Cl* > 0, with equality iff 
V = £ 2 = 0. However, since D, < 1 is bounded due to Eq. (jSHJ), we deduce that 
lim^+oo Q* = and, hence, that 

lim Q = 1, lim V = 0, lim £ 2 = 0. (42) 

T— »+0O T— > + 00 T— > + 00 

Thus, Eq. (I40J) implies that lim f ^ +00 g = (37 — 2)/2, from which it follows via Eqs. 
(jnHJ)-((SSJ) that for a sufficiently large f, there exists an e > such that dhiN a /df < —e 
and therefore that lim f ^ +00 iV a = 0. Consequently, a Bianchi type IX cosmology sourced 
by an ultra-stiff fluid isotropizes in the same way as the other Bianchi types as it 
approaches the big crunch at f — >• 00. In particular, this implies that there is no chaotic 
(oscillatory) behaviour in the vicinity of the singularity. 

Since we have now covered all possible Bianchi types, we may summarize the 
above analysis in the form of a cosmic no hair theorem: all initially contracting, 
spatially homogeneous, orthogonal Bianchi type I-VIII cosmologies and all Bianchi type 
IX universes that are sourced by an ultra- stiff fluid with an equation of state such that 
(7 — 2) is positive- definite collapse into an isotropic singularity, where the sink is the 
spatially flat and isotropic FRW universe. For the class of models where the equation 
of state asymptotes to a constant value 7 > 2 on the approach to the singularity, the 
FRW cosmology is a self-similar power-law solution, where the scale factor varies as 
a oc (-tf/^. 

In the following Section, we employ this theorem to determine the nature of a 
collapsing Bianchi type IX universe with a matter source consisting of a minimally 
coupled scalar field self-interacting through a negative- definite exponential potential. 



3. Isotropization of Collapsing Bianchi type IX Scalar Field Cosmology 

In this Section, we consider an action of the form 

S = Jd 4 x^g[R- OvV) 2 - V(<p)] , V = Vbexp(-Ap), (43) 

where V represents the potential of the scalar field, ip, and Vo < and A > are 
constants. We assume the field is orthogonal, i.e., that it is constant on the surfaces of 
homogeneity, cp = (p(t), so that its energy density is given by p = \(p 2 + V. The effective 
equation of state of the field is defined by 

7 = ^ (44) 

and is therefore bounded such that 7 > 2. 

The presence of a negative potential energy implies that not all the terms on the 
right-hand side of the Friedmann equation (j3J) will be semi-positive definite. Thus, a 
given Bianchi model may undergo a recollapse and, consequently, the Hubble normalized 
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variables do not lead to a global compact phase space. Nonetheless, the theorem 
developed in the previous Section can be employed to determine the nature of collapsing, 
homogeneous scalar field cosmologies in the vicinity of the big crunch singularity. 

It is known that the collapsing Bianchi I model is stable to curvature and anisotropy 
perturbations near the singularity [14 . In this Section, we will consider the Bianchi type 
IX universe. In this case, variables parametrizing the evolution of the scalar field can 
be defined such that [T3] 



The evolution equations for these variables are then determined by the scalar field 
equation, p = —3H(p 2 , which itself follows as a consequence of energy- momentum 
conservation. We find that 

$* = (?- 2)M - (46) 



:i+4>*-^* 



0. (47) 



It also proves convenient to define a new bounded variable |15j : 

" s dtt- (48) 

where < d < 1. This evolves such that 

d* = (1 + q)Hd{d - 1) (49) 

and the big crunch singularity, H — > — 1, corresponds to D — > oo (d — > 1). 

The physical state of a Bianchi type IX scalar field cosmology is therefore given 
by the vector x = (d, H, £+, Ni, N2, N3, ^, 0) and the corresponding autonomous 
set of ODEs are Eqs. (EU)-© and Eqs. flUJ), (gTJ) and (jUJ), subject to the constraint 
equations (J28)l and (|S9*jl . The density parameter is Cl = ^ 2 — 6 2 and evolves according 
to Eq. (jUJ. 

We will assume that the universe is undergoing a collapse, H < 0, and that the 
field has positive energy density, Cl > 0. It is anticipated that this should represent the 
behaviour of generic solutions once the recollapse has set in. Since q > 0, Eq. (p£T|) 
implies that H* < and, consequently, Eq. (j4*T| implies that Q continues to grow 
monotonically. Thus, the energy density of the field remains positive, <ft 2 > 2\V\, and 
the sign of ^ is fixed on the approach to the big crunch. Moreover, since Cl > and 
q > 0, the variables V, tt} < 1 (due to Eq. and lim^oo H = -1. 

The cosmic no hair theorem of Section II will apply for this collapsing Bianchi type 
IX scalar field cosmology if it can be shown that 7 > 2 at the singularity. In principle, 
however, a scalar field may become dominated by its kinetic energy during the collapse, 
in which case 7^2 and 0^0. This implies that the shear may not necessarily 
become sub-dominant. In order to establish whether this possibility arises, we will first 
show that if lim^oo 7 = 2, the contracting Bianchi IX model asymptotes in the vicinity 
of the big crunch toward the spatially flat type I background, where Ni = N 2 = N 3 = 0. 
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Firstly, an argument similar to that presented after Eq. (|41j) implies that 
lim^oo V — 0. However, for an equilibrium point with 7 = 2, Eq. ()36|) implies that 
q = 2. It then follows from Eq. ()32j) that this equilibrium point will also have S± = 0. 
The definition (j21J) then requires that N 3 = N 2 or N\ = N 2 + N 3 . If the latter condition 
applies, the requirements that V = S + = then imply that N 2 = (assuming without 
loss of generality that N 3 > N 2 > iV x ). On the other hand, if N x = N 3 ^ 0, Eqs. (0 
and (J33j) imply that E + = —1/2 and £_ = — \^3/2. In other words, the equilibrium 
point would correspond to a vacuum model, where ^ = O = due to the constraint 
()39|) . This is inconsistent with the property that Cl* > as f — > +00. 

Imposing the alternative condition N 2 = N 3 reduces the constraint S + = to 
N^Nx - iV 2 ) = 0. Hence, either jVi = N 2 = N 3 or JVi = 0. However, if N 2 = N 3 ^ 0, 
Eqs. (I34|) and (|35|l imply that S_ = and S + = 1, and once more this would correspond 
to a vacuum model, which is a contradiction. 

Hence, the only equilibrium point (with H < 0) that corresponds to a non-vacuum 
cosmology with 7 = 2 lies in the Bianchi type I invariant set and is given by 

d=l, H=-l, N 1 = N 2 = N 3 = 0, 

q = 2, 6 = 0, S^ + S 2 +^ 2 = l. (50) 

A standard perturbation analysis reveals that the eigenvalues associated with this point 
are 

-3, -4, 0, 0, -2-4S+, -4, 
-2 + 2S+±2 v / 3S_, -3-^-f. (51) 

The two vanishing eigenvalues imply that this is a two-dimensional set of equilibrium 
points. It follows that the stability of this set of points is sensitive to the kinetic energy 
of the scalar field. If ^ < — V~6/\, it represents a saddle, otherwise it may be a source 
if S + > 1/2 and S_ is sufficiently small. 

As there are no other equilibrium points, this implies in effect that if the scalar 
field is rolling to large negative values down its potential at a sufficiently fast rate, 
^ < — V6/X, its equation of state will satisfy lim f ^ 00 7 > 2. We may therefore apply 
the cosmic no hair theorem of the previous Section in this region of parameter space. 
In other words, the stable attractor solution is the spatially flat and isotropic FRW 
cosmology corresponding to the equilibrium point ^ = — X/V6, = [(A 2 /6) — l] 1 / 2 and 
q = (A 2 /2) — 1. This 'scaling' solution exists only for A > \^6. In conclusion, therefore, 
a contracting Bianchi IX cosmology sourced by an orthogonal scalar field rolling down a 
negative exponential potential with ^ < — V6/X and A > collapses into an isotropic 
singularity, where the sink is the self-similar, spatially flat and isotropic FRW universe 
with scale factor a oc (— t) 2 ^ 2 . 

In the next Section we consider such a scalar field model generated from a 
compactified higher-dimensional cosmology inspired by M*-theory. 
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4. Collapsing M*-Cosmology 



The M*-theory of Hull is a version of M-theory in an eleven-dimensional spacetime 
of signature (9+2). Compactification of this theory on a timelike circle gives rise to a 
string theory - termed the type IIA*-theory - with a supergravity limit that is similar 
to the conventional type IIA string but with a Ramond-Ramond (RR) sector containing 
form-fields whose kinetic terms have the 'wrong' sign. In this sense, M*-theory may be 
interpreted as the strongly-coupled limit of the type IIA* theory. Its low-energy limit is 
a supergravity theory with a bosonic sector given by 



5, 



d 



ii 



X \ 



(52) 



where G4 represents the field strength of the three-form potential C3. The sign of the 
kinetic term for G4 is determined by supersymmetry. 

Leaving aside issues associated with the conceptual problems of extra timelike 
dimensions (see jS] for a full discussion of such questions within the context of M*- 
theory), we may consider the Kaluza-Klein compactification of the action (JH2|) on a 
timelike circle S 1 . For the case where only the RR four- form field strength is non- 
trivial, the ten-dimensional action for the truncated type IIA*-theory is given by 



IIA* 



d x 



3-010 



R s + (W 



10 



G\ 
48 



(53) 



where the ten-dimensional dilaton field, 



no, 



is related to the radius of the eleventh 



dimension, e m , by 

0) 
9ab = 



ho = 3rn 
£l 2 gAB, 



and we have performed a conformal transformation: 
n 2 = e ri1 (54) 



(s) 

to the string- frame metric, g A g. 

We will now dimensionally reduce the theory ()53|) to four dimensions on a six-torus, 
T 6 , where the only dynamical degree of freedom in the internal dimensions is taken to 
be the breathing mode, (3. In other words, we assume that the string-frame metric (|54j) 
is given by 



dst 



9. 



${x)dxPd3? + e w{x) dy 2 , 



(55) 



where dy 2 = 5ijdy l dy J . We will also Hodge dualize the four- form field strength in ten 
dimensions to a six-form, Fq, and assume that the only non-zero components of this 
six-form live on the internal dimensions, i.e., we assume an ansatz Fq = mEQ, where 
Eq is the volume- form of T 6 and m is a constant. Thus, the effective four-dimensional 
effective action takes the form 



S = J d 4 x^f s \e-^ (R s + (V0 4 ) 2 - 6 (V(3f 



1 2 - 
—m e 
2 



•6/0 



(56) 



where 04 = <pio — 6/3 represents the four-dimensional dilaton field. 

Eq. may be expressed in the Einstein- Hilbert form by performing the conformal 
transformation 



g (e) 



Pi n {s) 
iJ (11/ 



(57) 
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and field redefinitions 4 = 4 / \/2 and (3 = \/Q[3. It follows that the conformally 
transformed action is given by 



S = Jd 4 x y^f e \R e - (V0 4 ) 2 - (V^) 2 + 
Finally, defining a new pair of scalar fields: 



(58) 



X 



f = (v/604 + VSP) (59) 



'14 

implies that the action (p)%|) is equivalent to 

S = Jd 4 x ^f e [R e - (Vx) 2 - (VO 2 + m 2 e-^] . (60) 

Thus, the potential for the x-field is negative and sufficiently steep (A = y/lA) 
for the results of Section III to apply The scalar degree of freedom £ will behave as 
a stiff perfect fluid, but the presence of this field does not affect the conclusions of 
the previous section, since it is decoupled from the other matter degree of freedom. 
It therefore becomes dynamically negligible near to the singularity in the case where 
the collapsing Bianchi type IX cosmology asymptotes to the FRW solution. In this 
case, the attractor is the power law solution a e oc (— t e ) 1 ^ 7 , \ — (2/\/l4) ln(— t e ), £ = 0. 
Transforming back into the four- dimensional string frame via Eq. (J57)) then implies that 
a s oc (— t s )~ 1//5 , e^ 4 oc (— t s )~ A ^ and oc (— t s ) 1 ^. The ten-dimensional string-frame 
metric ()55j) is therefore given by 

ds 2 s = T-^dxl + T l J 2 [-dT 2 s + dy 2 ] , (61) 

where we have defined T s = (— t s ) 4//5 and rescaled the coordinates where appropriate. 

The scale factors in the eleven-dimensional frame, which we may denote by 
{a, e@, e rn } are related to the corresponding ten-dimensional quantities by Eq. (|H1|) 
and take the form a = ae~ Tl1 ^ 2 and (3 = ft — rn/2, where the eleven-dimensional proper 
time is given by rj = J dt s exp(— rn/2). Hence, the eleven-dimensional line element can 
be expressed as 

ds 2 n = ff 2lz dx\ + r/ 1/3 [-dff 2 - dt 2 u + dy 2 ] , (62) 

where fj = (—i]) 6 ^ 7 , the coordinate of the eleventh (timelike) dimension is denoted by 
in and, for simplicity, we have absorbed any constants of proportionality by a suitable 
rescaling of the spacetime coordinates. 

A background of the form ()61)1 can be interpreted as the cosmological analogue of 
a domain wall spacetime . More specifically, we may write the metric in the form 
ds n = F~ 2 / 3 dxl + F 1//3 [— dr) 2 — dt\ x + dy 2 ], where {fj, tu,y} represent the transverse 
coordinates. Since F is a linear function of fj, it may be viewed as a harmonic function 
on the transverse space, thereby allowing Eq. (Jo^j) to be interpreted as a membrane- type 
solution with a three-dimensional Euclidean world- volume [TTH IT7] . In particular, it is 
known that solutions to M*-theory of this form preserve 16 of the 32 supersymmetries 
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of the theory ^7j. It is interesting that such a supersymmetric, higher-dimensional 
cosmological background is selected from a dynamical point of view when analyzed in 
terms of an effective four-dimensional, anisotropic cosmology. 

5. Conclusion 

To summarize, it has been shown that any orthogonal, collapsing Bianchi cosmology 
with a matter sector comprised of an ultra-stiff fluid with an arbitrary and varying 
equation of state (subject to the condition that 7 > 2 at all times) approaches spatial 
isotropy and flatness in the neighbourhood of the big crunch singularity. Of particular 
interest is the behaviour of the general Bianchi type IX universe. Since the singularity 
is isotropic, such models do not exhibit chaotic-type (Mixmaster) behaviour near to the 
big crunch, as is the case when 7 < 2 (see, e.g., ^H] and references therein). A related 
effect has been found in a class of braneworld models with a constant equation of state 
parameter 7 > 1, where the effective Friedmann equation depends quadratically on the 
energy density [TUj . This can be understood since the fluid in such a model behaves at 
small spatial volumes as if it was ultra-stiff [2*U] . 

The only spatially homogeneous cosmology that we have not considered is the 
Kantowski-Sachs model, where the G3 group of isometries acts multiply-transitively. 
The line element for this cosmology may be written in the form 

ds 2 = -dt 2 + D\dx 2 + D\dVL\ 

D 1 (t) = eftW-^W D2 (t) = e Mt)+Mt) (63) 

where dQf, is the metric on the two-sphere. Since this spacetime has positive spatial 
curvature, an extension of our cosmic no hair theorem to include this case would require 
a set of variables that compactifies the phase space in an analogous way to the variables 
(|26|) and (|2T|) for the Bianchi type IX model. However, for 7 > 2, a suitably normalized 
variable for the energy density that increases monotonically during a collapsing phase 
has yet to be identified. 

Nonetheless, we may gain insight into the asymptotic nature of an ultra-stiff 
Kantowski-Sachs universe by assuming 7 = constant. Following |21j, we define 
B12 = D\\- The Friedmann equation then takes the form 

3H 2 = p+±o»-B* (64) 

where H = O and a + = 3/3+ . It follows that D = [9H 2 + 3B 2 ] 1 / 2 is a dominant variable 
(assuming p > 0), and this implies that compact variables Qo = 3H/D, Q + = o~ + /D 
and Q = 3p/D 2 = 1 — Q 2 + may be introduced, together with a new time variable 
r = JdtD/3. The curvature variable is K = B 2 /(3H 2 ) = (1 - Ql)/Ql from which we 
deduce that < {Qo,Q + ,Q} < 1. 

It can then be shown j2I] that the field equations reduce to the set of evolution 
equations 

dQ 



dr 



l-Qt)F (65) 
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dQ+ 



i + (Qo - Q+f + QoQ+i 71 



(66) 



dr 



where 



F = y-l-QoQ + + ^(2- 7 )Q+ 



(67) 



Eq. then implies that all equilibrium points are located at Ql = 1 or F = 0. 
Let us first assume Ql ^ 1 and F = 0. It follows from Eqs. (|67)|) and (|H7|) that 
Q + = ±(37 — 2)/(4 — 37), but these points are unphysical for 7 > 2. Since we are 
interested in collapsing universes, we therefore consider the case Qo = —1- In this case, 
there are three equilibrium points, where Q + = 0, ±1. A stability analysis indicates 
that the points Q + = ±1 are saddles and the point (Q ,Q + ) = (—1,0) is a sink. This 
latter point corresponds to the isotropic, spatially flat FRW cosmology. Hence, in the 
case of a constant, ultra-stiff equation of state, the unique attractor for a collapsing 
Kantowski-Sachs cosmology is the isotropic and spatially flat universe. This provides 
strong support for the conjecture that this background represents the attractor for a 
general ultra-stiff equation of state. It also represents strong motivation for searching 
for an appropriate set of variables along the lines of Eqs. (|2T)|) - (j2*7|) . 

Our no hair results are also of relevance to the nature of spacetime singularities in 
more general inhomogeneous backgrounds, since according to the Belinskii, Khalatnikov 
and Lifshitz (BKL) conjecture [22], each spatial point of an inhomogeneous universe 
evolves effectively as a Bianchi model on the approach to the singularity. In principle, 
the results derived above can be employed to gain further insight into the nature of 
inhomogeneous cosmologies 0123 121] • Indeed, similar issues to those considered in this 
paper have been addressed by Coley and Lim [23; who performed an asymptotic analysis 
of stiff and ultra-stiff Abelian G 2 and general Go spatially inhomogeneous cosmologies, 
and found that the spatially flat FRW solution is locally stable in the past for 7 > 2. 
Our results differ in that it was assumed that 7 is a constant in [21], whereas we have 
assumed an arbitrary (differentiable) equation of state 7 = 'j(p) > 2. Furthermore, our 
no hair result for the spatially homogeneous Bianchi models is a global result, in the 
sense that the flat FRW solution is the unique attractor at the big crunch. 

Finally, the work of [21] is restricted to that of a perfect fluid matter source. 
We extended our analysis to the case of a collapsing, spatially homogeneous Bianchi 
type IX model containing a scalar field that self-interacts through a negative potential. 
In this case, the effective equation of state varies but satisfies 7 > 2, with equality 
corresponding to a vanishing potential. Our cosmic no hair theorem implies that the 
problem of demonstrating that such models isotropize as they collapse has been reduced 
to verifying that the potential remains dynamically significant near to the singularity. 
We have shown that this is possible, for example, in the case of a steep, negative 
exponential potential. Exponential potentials are of particular interest since they are 
common in cosmological models inspired by string/M-theory. One such theory that 
leads naturally to negative potentials is M*-theory, where the potential is generated 
by non-trivial flux of the form fields |5\. We considered a truncated, dimensionally 
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reduced effective action for this theory and found that when the isotropic attractor 
for a collapsing Bianchi type IX cosmology is oxidized back to eleven dimensions, the 
metric is a supersymmetric solution of the higher-dimensional theory. This is of interest 
since supersymmetric backgrounds are non-perturbatively exact solutions and therefore 
provide a valuable framework for investigating the cosmic dynamics in the high curvature 
regime near to the singularity. We expect that other attractor solutions for more general 
dimensional reduction schemes will also exhibit similar supersymmetric properties when 

interpreted in a ten- or eleven-dimensional context. 
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